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Abstract We review our recent results of Alfven wave-driven winds. First, we present 
the result of a self-consistent ID MHD simulations for solar winds from the photo- 
sphere to interplanetary region. Here, we emphasize the importance of the reflection 
of Alfven waves in the density stratified corona and solar winds. We also introduce 
the recent HINODE observation that might detect the reflection signature of trans- 
verse (Alfvenic ) waves by Fujimura & Tsuneta (2009). Then, we show the results of 
Alfven wave-driven winds from red giant stars. We explain the change of the atmo- 
sphere properties from steady coronal winds to intermittent chromospheric winds and 
discuss how the wave reflection is affected by the decrease of the surface gravity with 
stellar evolution. We also discuss similarities and differences of accretion disk winds by 
MHD turbulence. 

Keywords accretion disks ■ magnetohydrodynamics ■ solar wind • stellar winds • 
turbulence • waves 



1 Introduction 

The main difficulty of the coronal heating and solar wind acceleration is how to lift 
up nonthermal (e.g. magnetic) energy from the photosphere to upper regions because 
the hotter corona cannot stably exist above the cooler photosphere by upward thermal 
heat flux. In addition, it is not well understood how to let the energy dissipate at 
appropriate positions. The origin of energy that heats up the corona and accelerates 
the solar winds is in the surface convective layer. The solar atmosphere is filled with 
complicated structure of magn etic field which is supposed to be amplified by dynamo 
mechanism in the interior (e.g.. lBrun et al.ll2004 ). The turbulent motions of the surface 
convection drive various modes of upward propagating waves. The turbulent motions 
may also trigger magnetic reconnections which result in flares and flare-like events, 

Takeru K. Suzuki 

Puro-cho, Chikusa, Nagoya, 464-8601, Japan 
Tel.; -1-81-52-788-6196 
Fax; +81-52-788-6196 
E-mail; stakeru@nagoya-u.jp 



2 



which drive various mode of waves at locations above the photosphere l|Sturrocklll999l l. 
In this paper, we focus on the roles of the waves in heating and accelerating the solar 
winds. 

The compressive waves that are excited at the photosphere cannot travel to suffi- 
ciently upper regions because they easily steepen the wave fronts to form shocks after 
the amplification of the amplitude in the density decreasing atmosphere. For instance, 
the acoustic waves mostly dissipate before reaching the corona, and therefore they 
cannot contribute to the heating of th e corona and the acceleration of the solar wind 
( Stein fc Schwartz II 19721 : ISuzukill2002h . Fast mode waves in low /3 (magnetically dom- 
inated) condition suffer refraction so that they hardly reach the coronal height if the 
fast mode speed varies in ho rizontal direction due to complicated magnetic structure 
(|Khomenko fc Colladoslbooel ). In general, the compressive waves that are excited from 
the photosphere cannot give a significant contribution to the coronal heating and solar 
wind acceleration. 

On the other hand, Alfven waves can travel a longer distance owing to the incom- 
pressive characters. A fraction of the Alfven waves excited at the photospheric level 
is supposed to propagate to the corona and the solar wind acceleration region so that 
they play a more dominant role in the heating and acceleration of the solar wind than 
the compressive waves. Alfvenic oscillations are also detected in various regions on the 
Sun by the recent obs ervations usi ng the HINODE satellit e, e.g., in a solar prominence 
(Okamoto ct al. '2007*), in spicules dPe Pontieu et aTll2007l l. and in a chromospheric jet 
(Nishizuka ot al. 2008J. 

To understand the dissipation mechanism of the Alfven waves is a key to understand 
the heating and acceleration of the solar wind because this corresponds to the transfer of 
the energy flux of the Alfven waves to the ambient plasma. The Alfven speed in general 
largely changes wit h height because th e density and magnetic field strength decrease in 
a different manner. IShiota et al.l ( 2010l ) have estimated the distribution of Alfven speeds 
in a polar region from HINODE observation and found that the Alfven speeds ac- 
tually vary in both vertical and horizontal directions. Under these circumstances, 
Alfven waves suffer reflection through the upward propagation as a result of the de- 
formation of the wave shape (|An et al.|[l99ol : I Moore et al.|[l99ll : iHollweg fc Isenberd 
Interactions between the outgoing Alfven waves and the i ncoming com p onent 
possibly lead to the damping through turbu lent-like cascade (e.g.. iDmitruk et al.ll2003l : 
ICranmer. et ^l2007l : IVerdini fc VeUiiboOTl ') . This process might further generate high- 
frequency ioncyclotron waves, which are widely discussed especially on the preferential 
heating of the perpendicular components with resp ect to magnetic field of minor heavy 
ions fe.g.. E\i fc Marschlfl997l : iKohl et"aHll998l : [Cr"anmer. et al.lll999i '). The variation of 
the A lfven speeds also anticipate phase mixing of Alfven waves along neighboring fiel d 
lines (|Hevvaerts fc PriestI [l983 : ISakurai fc Graniklll984l: IVoitenko fc Goosseiii 120051 ) . 



which might further generate fast mode waves (jNakariakov et al.lll997l ). The density 



stratification and the variation of the Alfven speed also enhances the parametric decay 
instability of outgoing Alfven waves (Suzuki fc Inutsuka 2006; SI06 hereafter). The 
parametric decay excites compressive waves in addition to incoming Alfven waves. In 
contrast to those excited at the photosphere, these compressive waves are generated in 
the upper regions and directly contribute to the heating of the corona and solar winds 
by the shock dissipation. 

Although the final wave dissipation mechanism and heating process are still not 
well figured out, understanding the wave reflection in stratified atmosphere with vary- 
ing Alfven speed is primarily important. Alfven waves are refiected most effectively 
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in the chromosphere and transition region because of the rapid change of the density. 
While most of the Alfven waves from the surface cannot penetrate into the corona, 
the remaining waves sufficiently heat and accelerate the solar wind. Our previous work 
(Suzuki & Inutsuka 2005; SI05 hereafter) showed that the transmitted fraction into the 
corona is an order of 10% and that this is enough for the coronal heating and solar 
wind acceleration in open coronal holes. Although turbulent-like cascade of Alfven (ic) 
waves is not treated in this work because we adopted one dimensional (ID) magneto- 
hydrodynamical (MIfD) simulations, this is the first attempt that dynamically treats 
wave reflection from the photosphere to solar wind regions simultaneously with heat- 
ing of the gas by the wave damping. In this review paper, we flrstly summarize our 
dynamical simulations of solar winds, focusing on the reflection of Alfven waves under 
various conditions. Actua lly, the recent HINODE ob servation also detected a signature 
of reflected Alfven waves ( Fuiimura fc TsunetallioogI ) . which will be introduced in more 
detail later. 

The concept of Alfven wave-driven winds is not limited to the S un. It is ex pected 
that the simil ar mechanism o perates in winds from other stars fe.g. IVeliilll993l \ such 
as protosta rs (ICran mei| l200^ ). low- and intermediate-mass main sequence stars, red 
giant stars (^Suzuk?200r, S07 hereafter), and proto- neutron stars with strong magnetic 
held (Suzuki & Nagataki 2005). Among these objects, we show Alfven wave driven red 
giant winds in this paper based on our work (S07). In red giant stars, the stratification 
properties (e.g. decrease of density) of the atmosphere are different from those in the 
Sun because of the smaller surface gravity. This affects the propagation and reflection 
of the Alfven waves. Here we discuss the change of the wave reflection with stellar 
evolution. 

In the last part of this paper, we further extend to winds from accretion disks. To 
date, i t has been widely discussed that accretion disk wind s can be driven by centrifugal 
force l|Blandford fc PavnjflQsi : iKudoh fc Shibatalll998h . Global magnetic field plays 
a key role in this process; if there is poloidal field that is sufficiently tilted with respect 
to an accretion disk, the gas can flow out along with the held lines. On the other hand, 
it is expected that small-scale MHD turbulence also potentially drives disk winds, 
similarly to the solar winds, although such mechanism has not been well studied so 
far. It is now widely accepted that the outward transport of disk angular momentum 
and inward mass accretion are realized by the effective viscosity owing to turbulence. 
Magnetorotational instability (MRI) is now known to be a very effective mechanism 
that drives turbulence if weak magnetic fields exist initially (Balbus fc Hawley 1991). 
The excited MHD turbulence is supposed to accelerate the disk material to upper 
regions as well as transport the angular momentum outwardly. In Suzuki fc Inutsuka 
(2009), we studied the disk winds driven by MHD turbulence by local disk simulations. 
In this paper, we introduce the results with emphasizing the similarity to and difference 
from the Alfven wave-driven solar and stellar winds. 



2 Solar Winds 

2.1 Simulation Set-up 

We simulate the heating and acceleration of solar winds in ID super-radially expansing 
flux tubes from the photosphere to sufficiently distant locations (0.1 - 0.3 AU). Here 
we summarize the basic points of the simulation set-up. For the details, please see SI05 



Fig. 1 Set-up of solar wind simulations. The cartoon shows a schematic view of a super- 
radially open flux tube, which we set in our simulations. 



and SI06. We determine radial magnetic field strength, Br, from the conservation of 
magnetic flux : 

Brfr"^ = const., (1) 

where r is heliocentric distance, and / is a super-radial expansion factor. We give radial 
variation of / in advance. The structure of Brir) is determined by surface radial field, 
Bj-fi, and /, and is fixed during the simulations. In our simulations (SI05 & SI06), 
we consider the magnetic field of ~ few hundred Gauss for 5^,0 a^nd the total super- 
radial expansion factor of 75-450, which give the average magnetic field strength of 
open field regions, B^fl/ f , of an order of 1 G. We should remark that recent HINODE 
observations found that open magnetic flux tubes in polar regions ar e connected to 
ubiquitously distributed unipolar patches with more than kilo-Gauss ( Tsuneta el al.l 



boOSi : IShiota et al.ll2O10l ). These open flux tubes super-radially open with a factor of 
several hundred to a thousand. Although the obtained _Br,o a^nd / are slightly larger 
than those adopted in our sim ulations, the v alues of the average strength, -Br,o//, are 
similar. As shown in SI06 and ISuzukil 1 20061 ). B^.o/f is a more important parameter 



in determining the properties of the solar winds. 

We inject transverse perturbations of magnetic field with various spectra from the 
photosphere. Outgoing Alfven waves are excited by these fluctuations. The propagation 
and dissipation of the waves are dynamically treated by the following MHD equations 
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dVr dp 1 9 2.132 ^ , P"i 9 / 2j.x GMq 



Br 
47r ' 



+ ^^('^'/^c) + gR = 0, (5) 



where p, v, p, B are density, velocity, pressure, and magnetic field strength, respectively, 
and subscripts r and ± denote radial and tangential components; ^ and ^ denote 
Lagrangian and Eulerian derivatives, respectively; e — ^pri^ is specific energy and 
we assume the equation of state for ideal gas with a ratio of specific heat, 7 = 5/3; G 
and M© are the gravitational constant and the solar mass; Fc[= KpT ^^^^) is th ermal 



conductive flux by Coulomb collisions, where kq = 10 ^ i n c.g.s unit ( Braginskii 19 65|): 

)ri- Fossil 



f?R is radiative cooling. We use optically thin radiative loss (jLandini fc Monsignori- 
Il99dl in the corona and take into account opti cally thick effects in the chromosphere 
(|Anderson fc Athavlll989l : iMorivasu et al.ll2004l . SI06). 

We adopt the second-order MHD-Godunov-MOCCT scheme to update the physical 
quantities. Each cell boundary is treated as discontinuity, and for the time evolution 
we solve nonlinear Riemann shock tube problems with the magnetic pressure term by 
using the Rankin-Hugoniot relations. Therefore, heating is automatically calculated 
from the shock jump condition. A great advantage of our code is that no artificial 
viscosity is required even for strong MHD shocks; numerical diffusion is suppressed to 
the minimum level for adopted numerical resolution. 

We initially set static atmosphere with a temperature T = lO'^K to see whether the 
atmosphere is heated up to coronal temperature and accelerated to a transonic flow. 
At t = we start the inject of the transverse fluctuations from the photosphere and 
continue the simulations until the quasi-steady states are achieved. 



2.2 Results 

Before discussing the reflection of the Alfven waves in various cases, we explain how the 
coronal heating and the solar wind acceleration were accomplished in the typical case 
which we studied in SI05 and SI06 for the fast solar wind. We adopt B^fi ~ 161 G, the 
total / = 75, and the root-mean-squared (rms) surface amplitude, (dv) = 0.7 km s~^. 
Figure [2] plots the flnal structure of the simulated solar wind after the quasi-steady 
state is achieved in comparison with observations of fast solar winds. In the four panels 
(km s-l), r(K), mass density, p(g cm '^), and rms transverse amplitude, {dv±){km 
s~^) are plotted. As for the density, we compare our result with observed electron 
density, A^e, in the corona. When deriving A^e from p in the corona, we assume H and 
He are fully ionized, and Afe(cm~^) = 6 x 10^^p(g cm~^). 
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Fig. 2 Comparison of the simulation result with observations of fast solar wind. Outflow speed, 
^^(km s~^) (Top-left), temperature, r(K) (top-right), rms transverse amplitude, {dv±)(km 
s~^) (bottom- left), and density in logarithmic scale, log(p(g cm~^)) (bottom-right), are plot- 
ted. Observational data in the temperature panel are electron density, log(Afe (cm~^)) which 
is to be referred to the right axis. Dashed lines indicate the initial conditions and solid lines 
are the results at i = 2573 minutes. In the bottom panel, the initial value {{dvj_) = 0) 
dose not appear. Top-left: Green vertical error bars are proton outflow speeds in an inter- 
plume region by UVCS/SoHO (Teriaca ct al. 2003). Dark blue vertical error bars are proton 
outflo w speeds by the Doppler dimming technique using UVCS/SoHO data dZangrilli et al.l 
120021) . A dark blue open square with errors is velocity b y IPS measurements averaged in 
0.13 - 0.3A U of high-latitude regions llKoiima et al.ll2004l '). Light blue data are taken from 
iGrall et al.l (1996); crossed bars are IPS measurements by EISCAT, crossed bars with open 
circles are by VLBA measurements, and vertical er ror bars with open circles are data based 
on observation by SPARTAN 201-01 (Habbal et al. 1994). top-right: Pink circles are electron 
temperatures by CDS/SoHO l lFludra ct al. 1999). botto m-left: Blue circles are non-thermal 
broadening inferred from SUMER/SoHO measurements llBaneriee et al.lll99 8). Cross hatched 
region is an empirical co nstraint of non-thermal broadening based on UVCS/SoHO obser- 
vation (Es ser et al ] |1999D . Green error bars are tra nsverse velocity fluctuations derived from 
IPS measurements by EISCA Tl|C anal s et al. ir200^. bottom-right: Circl es and stars are ob - 
servations by SUMER/ SoHO l|W ilhclm ct al. 199^) and b y CDS/SoHO llTbriaca et al.ll2003) . 
respectively. Triangles llTeriaca et al.l 12003) and squares llLamv et al.lll997D are observations 
by LASCO/SoHO. 
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Figure [2] shows that the initially cool and static atmosphere is efltectively heated 
and accelerated by the dissipation of the Alfven waves. The sharp TR which divides 
the cool chromosphere with T ~ lO^K and the hot corona with T ~ lO^K is formed 
owing to a thermally unstable regio n around T ~ 10^ K in the radiative cooling function 



( Landini fc Monsignori-Fossilll99Cl ). The hot corona streams out as the transonic solar 



wind. The simulation naturally explains the observed trend quite well, (see SI05 and 
SI06 for more detailed discussions.) 

The reflection of Alfven waves are easily seen in r — t diagrams. Figure |3] presents 
contours of amplitude of Vr, p, v^, and Bj^/Br in Rq < r < 15J?© from t — 2570 min. 
to 2600 min. Blue (gray) shaded regions denote positive (negative) amplitude. Above 
the panels, we indicate the directions of the local 5 characteristics, two Alfven , two 
slow, and one entropy waves at the respective positions. In our simple ID geometry, Vr 
and p trace the slow modes which have longitudinal wave components, while vj_ and 
B_L trace the Alfven modes which are transverse (note that fast-mode and Alfven mode 
degenerate in the simple ID treatment, and then we simply call them Alfven waves). 

One can clearly see the Alfven waves in vj_ and Bj_/Br diagrams, which have the 
same slopes with the Alfven characteristics shown above. One can also find the incom- 
ing modes propagating from lower-right to upper-left as well as the outgoing modes 
generated from the surfac^. These incoming waves are generated by the reflection at 
the 'density mirrors' of the slow modes in addition to the reflection owing to the shape 
deformation (SI05 & SI06). At intersection points of the outgoing and incoming char- 
acteristics the non-linear wave-wave interactions take place, which play a role in the 
wave dissipation. 

The slow modes are seen in Vr and p diagrams. Although it might be difficult to 
distinguish, the most of the patterns are due to the outgoing slow mo defl which are gen- 



erated from the perturbatio ns of the Alfven wave pressure, BI/Stv (|Kudoh fc Shibatal 
1999l : lTsurutani et al.ll20o3 v These slow waves steepen eventually and lead to the shock 
dissipation. 

The processes discussed here are the combination of the direct mode conversion to 
the compressive waves and the parametric decay instability due to th ree-wave (outgo- 



ing Alfven , incoming Alfven , and outgoing slow waves) interactions ( Goldstei"nlll978l : 
iTerasawa etall 119861 of the Alfven waves. Although they are not generally efficient 
in the homogeneous background since they are the nonlinear mechanisms, the density 
gradient of the background plasma totally changes the situation. Owing to the gravity, 
the density rapidly decreases in the corona as r increases, which results in the amplifi- 
cation of the wave amplitude so that the waves easily become nonlinear. Furthermore, 
the Alfven speed varies a lot due to the change of the density even within one wave- 
length of Alfven waves with periods of minutes or longer. This leads to both variation 
of the wave pressu re in one waveleng th and partial reflection through the deformation 



of the wave shape ( Moore et al.lll991 '). The dissipation is greatly enhanced by the den- 
sity stratification, in comparison with the case of the homogeneous background. Thus, 
the low-frequency Alfven waves are effectively dissipated, which results in the heating 
and acceleration of the coronal plasma. 



It is instructive to note that tlie incoming Alfven waves have the positive correlation 
between and (dark-dark or light-light in the figures), while the outgoing modes have 
the negative correlation (blue-gray or gray-blue). 

^ The phase correlation of the longitudinal slow waves is opposite to that of the transverse 
Alfven waves. The outgoing slow modes have the positive correlation between amplitudes of 
Vr and p, {SvrSp > 0), while the incoming modes have the negative correlation [SvrSp < 0). 
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Fig. 3 r — t diagrams for Vr (upper-left), p (lower-left), v± (upper-right), and B±/Br (lower- 
right.) The horizontal axises cover from Rq to 15Rq, and the vertical axises cover from 
t = 2570 minutes to 2600 minutes, where Rq is the solar radius. Blue and gray shaded regions 
indicate positive and negative amplitudes which exceed certain thresholds. The thresholds are 
dvr = ±96km/s for Vr, dp/p = ±0.25 for p, v± = ±180km/s for v^, and B±/Br = ±0.16 for 
B± I Br , where dp and dvr are differences from the averaged p and Vr ■ Arrows on the top panels 
indicate characteristics of Alfven , slow MHD and entropy waves at the respective locations 
(see text). 
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Fig. 4 The top panels show the wave actions (5c) of the out-going (soUd) and incoming 
(dashed) Alfven waves. The bottom panels show density (solid) and temperature (dashed). 
The left panels are the results of the surface amplitude, {dv) = 1.4km s~^, the middle panels 
are the results of {dv) = 0.7km s~^, and the right panels are the results of {dv) = 0.4km s"'^. 



From now, we examine how the reflection properties are affected by changing in- 
put parameters. In order to quantitatively study the dissipation and reflection of the 
Alfven waves, we use an adiabat ic constant, Sc, of the outgoing Alfven wave derived 
from wave action (|jacaues Ill977l ): 

where vj^ = Br/ -^47rp is Alfven velocity and 

5«A,± = ^ =F B±/v/4^) (8) 
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is the amplitudes of the outgoing (for + sign) and incoming (for — sign) Alfven waves 
(Elsasser variables). 

First, we show the results with different values of surface fluctuations, (dv) = 
0.4, 0.7, 1.4 km s^^ (Figure^. In the upper panel, we compare Sc.+ and Sc,- of these 
three cases, whereas we use the absolute values for Sc,- to show in logarithmic scale. 
As shown in the bottom panels, the density is very sensitive to the input (dv). The 
density of the (dv) = 1.4 km s~^ case is 1000 times larger than the density of the 
(dv) = 0.4 km s""*^ case. Accordingly the mass fluxes of the solar winds differs about 
several hundred times. The differences of the densities and the mass fluxes between 
the two cases are much larger than the ratio of the input energy (oc dv'^) « 12. The 
reasons why the density sensitively depends on the input surface perturbation can be 
explained by the reflection and the nonlinear damping of Alfven waves. 

The difference of the wave reflection in different (dv) cases are illustrated in the 
upper panels of Figured In the smaller (dv) cases, the values of the outgoing (solid) 
and incoming (dashed) components are very similar especially in the chromospheric 
regions (r — Rq ~ O.OI-R0). This indicates that most of the outgoing Alfven waves 
generated from the surface is reflected back downward. Because the heating is smaller, 
the temperature is lower in the smaller (dv) cases. Then, the scale height becomes 
smaller and the density decreases rapidly. The Alfven speed changes more rapidly and 
the wave shape is largely deformed, which enhances the reflection. When the input 
wave energy decreases, a positive feedback operates; a smaller fraction of the energy 
can reach the coronal region. As a result, the density and mass flux of the solar wind 
becomes much smaller than the decreasing factor of the input energy. 

In addition to the effect of the wave reflection, the nonlinear dissipation of the 
Alfven waves also plays a role in the sensitive behavior of the density on the input 
wave energy (SI06). When the input wave energy becomes smaller, the density becomes 
smaller as explained above. Then, the nonlinearity of Alfven wave, 5ua,+/i'A) decreases 
not only because the amplitude, 5t;A,+ , is small but also because the Alfven speed, 
va(oc is larger. Therefore, the Alfven waves do not dissipate and the heating is 

reduced, which further decreases the density; this is another type of positive feedbacks, 
which also results in the sensitive dependence of the density on the input wave energy. 

Next, we investigate how the wave reflection is affected by the spectra of the input 
fluctuations. Here we study the three cases : The first case adopts the power spectrum, 
P{v) oc (pink noise), with respect to frequency, v; the second case adopts P{v) oc 
u'^ (white noise); the third case adopts sinusoidal perturbation with period of 3 minutes. 
In the first and second cases, the range of the period is set from 20 seconds to 30 
minutes. Figure [5] compares these three cases. The first and second cases show similar 
results, but in the chromospheric regions the first case suffers more reflection. This is 
because larger power is in longer wavelengths (lower frequency). If the wavelength is 
longer than the characteristic variation scale of Alfven speed, Alfven waves suffer more 
reflection. 

The third case (sinusoidal perturbation) shows considerably different structure from 
these two cases. The transition region is not so sharp, which is seen in the temper- 
ature structure (dashed line in the bottom panel), because the heating is localized 
near the wave crests which are formed as a result of steepening of nonlinearly excited 
compressive waves. Then, the heating becomes more sporadic and responding to this 
intermittent heating the transition region moves up and down. This is a clear contrast 
to the first and second cases in which the sharp transition regions work as walls against 



11 



10' 



I 

to 

CM 

u 

CD 



CO 



0^ r 



10 

10" 

10" 

10" 

^10" 

'^10" 

<10" 
cn 

^10" 
10" 
10" 
10" 



I I llllll| — I I llllll| — I I llllll| — I I llllll| H-l I llllll| — I I llllll| — I I llllllj — I I llllllj H-l I llllll| — I I llllll| — I I llllll| — I I llllll| 




-. 10° 



10^S~ 



- 10^ 



0.01 0.1 1 10 0.01 0.1 1 10 0.01 0.1 1 10 

(r-R0)/R0 (r-Ro)/Ro (r-R0)/R0 

Fig. 5 Same as Figure 2] but for the cases with different spectra of surface perturbation. 
The left panels are the results of P(t') oc (pink noise), the middle panels are the results 
of P{u) oc (white noise), and the right panels are the results of sinusoidal waves with 3 
minutes. 



the Alfven waves because the Alfven speeds change drastically. On the other hand, in 
the sinusoidal case Alfven waves become more transmittant owing to the fluctuating 
transition region. As shown in the top panel, the reflected fraction (the ratio of the 
incoming component to the outgoing one in the chromosphere) is smaller in this case. 



2.3 Discussions 

2.3.1 Recent Observation of Waves 

Wave activities are observed in various portions of the solar atmosphere. Propagat- 
ing slow MHD waves are often observed in the corona as the fluctuations of intensity 
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( Ofman et aHll999l : ISakurai et al]|2002l : iBaneriee et~aH I2OO9I I . iTomczvk et~aH (|2007l ') 



detected Alfven waves propagating along with magnetic field lines in the corona. Al- 
though the obtained wave amplitude is smaller than t he value required for th e coronal 
heating, this may be due to the line-of-sight effect I Mcintosh et alj l2008h . Indeed, 



nonthermal broadening of spectral lines, w hich are supposed to be closely linked to 
Alfven waves, gives much larger amplitudes I Singh et al.ll20o3 : IBaneriee et al.ll2009h 



Alfven (ic) waves are detected at the photospheric level in quiet regions (jUhichl 

Il996l l by analyzing the phase correlation of magnetic field and veloc ity perturbations . 
Similar attempts have b een carried out in sunspot regions as w ell ( Lites et al] Il998l : 



iBellot Rubio et all I2OO0I ) , Recently, iFuiimura &: Tsunetal ( 20091 ) carried out very de- 



tailed analysis of the observ ed fluctuati ons by the spectro-polarimeter ( SP) of the So- 
lar Optical Telescope (SO T; iTsuneta et al. 2008; Sucmatsu ct al. 2008; I chimoto et al.l 
I2OO8I : IShimizu et al.lliooj 'l aboard the HINODE satellite. They obtained the fluctua- 
tions of magnetic field, SB, the fiuctuations of velocity field, Sv, and the fiuctuations 
of intensity, SI, which reflects density perturbation, Sp of pores, and magnetic concen- 
trations in plages. The periods of the oscillations are distributed from 3 to 6 minutes 
in the pores and from 4 to 9 minutes in the plages. From the phase correlation of 
SB, Sv, and SI they investigated the modes and directions of the waves to that are 
attributed the observed fluctuations. One of the most important results is that the 
phase differences between Sv and SB are nearly -90 degrees 0. The obtained phase 
correlation is consistent with standing kink-mode (transverse) or sausage-mode (longi- 
tudinal) waves, which are, roughly speaking, surface-mode counterparts of Alfven and 
slow MHD waves. Namely, waves that propagate to positive and negative directions 
are almost equally exist at the photospheric level. Although it is quite difficult to es- 
timate the fraction of sausage (compressive) mode and kink (nearly incompressive) 
mode because the transformation from SI to Sp is not straightforward, it is expected 
that regions with smaller SI are dominated by kink-mode oscillations. Based on this 
consideration, they picked up a particular pore region with relatively small SI. In this 
region, the phase difference between SB and Sv is -96 degrees, which indicates that the 
upward propagating flux is slightly larger than the downward propagating component. 
Assuming the observed fluctuations attribute to kink-mode waves, they obtained the 
net Poynting flux that propagates upward as 2.7 x 10® erg cm^^s^^. 

Interestingly enough, our simulations also show very effective reflection of Alfven waves 
below the transition region. In the reference case for the fast solar wind (Figure [5]), 
85% of the initial upward Poynting flux of the Alfven waves from the photosphere is 
reflected back downward before reaching the corona. This implies that one can observe 
the comparable amount of downward flux to the upward component. The net leakage of 
the upgoing flux to the corona in this simulation is 15% of the input at the photosphere 
(~ 5 X 10^ erg cm^'^s"^). These f eatures seem quite consisten t, at least in a qualita- 
tive sense, with the observation bv lFuiimura fc Tsunetal (|2009h . The specific values of 
the leaking and reflected fluxes depend on strength and configuration of magnetic flux 
tubes and amplitudes of surface fluctuations. 



^ IFuiimura &c Tsunetal ll 20091) picked up tfie regions witli positive magnetic polarity (mag- 
netic fields toward observers), whicli is tiie reason why there are no data with the phase shift 
of -1-90 degrees. Data with negative polarity actually show the phase shift peaked around -|-90 
degrees (Fujimura & Tsuneta, private communication). 
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2.3.2 Data Driven Simulations 

The observation data are also used in numerical simulations. Recently, various groups 
are perfo rming so-called data-driven simulations by ado pting magnetic field data on the 
Sun fe.g. [Manchester et ahlfioOSl : iKataoka et al.ll2009l ). While these works are mainly 
aiming at global heliospheric phenomena such as propagation of coronal mass ejections, 
observed data by HINODE a re also applied for MHD simulations of surface activities. 
iMatsumoto &i Shibatal ( 20091 ) have carried out numerical simulations in open field re- 



gions by using observed transverse motions on the Sun. They obtained horizontal mo- 
tions of granules by local correlation tracking|f|. These transverse motions are expected 
to excite Alfvenic waves that propagate along the vertical magnetic fields. They ob- 
served 14 different regions and derive the power spectrum of the transverse oscillations. 
Using the obtained spectrum as the input from the photosphere, t hey performed ID 
MHD simulations in open field regions up to the corona following iKudoh fc Shibatal 
( 19991 ). They found that, compared to artificial inputs of white noise or pink noise, more 



fraction of the Alfven wave energy transmits through the transition region. They inter- 
pret that this is because more energy is resonantly trapped between the photosphere 
and the transition region and eventually leaks into the corona. 



2.3.3 Slow Solar Winds 



In this paper, we have introdu ced our simula tion results for the fast solar winds from 
polar coronal holes. In SI06 and lSuzukil ( 200(i ) . we tried to explain observed data of slow 



solar winds in the same framework of the nonlinear dissipation of the Alfven waves. 
We have shown that the slow winds from fiux tubes with smaller Bj.fi/ f and slightly 
larger {dvj_) well explain the observed properties of the slow streams from mid- to 
low-latitude regions; the same physical mechanism operates but the environments (e.g. 
geometries of fiux tubes) make the differences between the fast and slow solar winds. 
This interpretation is q uite consistent with the observed data by Interplanetary scin- 
tillation measurements ( Koiima et ahllioosl ). 



HINODE o bservations detected mass outfiows from open regions i n the vicinities of 
active regions (jSakao et al.ll2007l : llmada et aPboOTl : iHara et al.ir2008h . These outfiows 



might become slow solar win ds and possibly giv e the significant contribution to the 



total mass loss from the Sun i Sakao et al.ll2007l '). What is puzzl i ng is that the accel 



eration seems to take place at a very low altitude. iHarra et al.l l|200a) measured the 
Doppler velocities of the outfiows and concluded that the outfiow speeds reach ~ 100 
km at the height of ~ O.IJ?© from the surface. Such rapid acceleration is very 
di fferent from those infe rred in the classical slow winds. For instance, the observations 
bv ISheelev et al. 1 lll99'itl shows that the slow winds associated with the streamer belt 



reach ~ 100 km s ^ above 2-3 solar radii. 

It is difficult to explain the observed rapid acceleration by our simulations of the 
nonlinear Alfven waves. If we increase the input energy (~ (dv'^)), density, rather than 
outflow speeds, increases owing to the larger heating. Probably other calculations or 
simulations will give similar tendency. To explain the outflow from the near active 
regions, pure momentum inputs are necessary. 



* To do so HINODE/SOT is the best telescope, because horizontal oscillations are largely 
affected by the atmospheric seeing. 
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3 Evolution to Red Giants 



When the Sun evolve to a red giant star, the properties of the stellar wind also change. 
When the stellar atmosphere becomes cool enough to form dusts during the asymptotic 
giant branch ( AGB), the ste llar wind is mainly driven by the radiation pressure on the 
dust particles I Bowenlll988l ). Before the AGB phase, namely from the main sequence 
to (early) red giant branch (RGB), the driving mechanism of the stellar winds is re- 
garded to be similar to the solar wind; the origin of the driver is the kinetic energy 
of the surface convection, and MHD waves generated by the turbulent motion of the 
surface convection ac celerates the stel lar winds (e.g. Hartmann fc MacGregor I ll980l : 
Ijatenco-Pereira fc Qpher.lQS&.Charbonneau fc MacGregor iggSl. However, the prop- 
erties of the stellar winds change with stellar evolution primarily because the surface 
gravity decreases. 

Based on these considerations, we studied the evolution of the stellar winds with 
the stellar evolution from the main sequence to the RGB (S07). We extend the solar 
wind simulation introduced in the previous section to simulate the red giant winds. We 
change the stellar radius, which controls the surface gravity, and effective temperature 
(or sound speed). We input the surface fluctuations that excite outgoing waves from 
the photosphere in the same manner as in the solar wind simulation. The amplitudes 
and spectra of the perturba tions are estimated from the scaling relation of surface 
convective flux and period faenzini et al.l Il977l . see also Brun & Palacios 2009 for 
recent simulations). We simulated the stellar winds from a IMq star with different 
surface gravities, logg = 4.4, 3.4, 2.4, and 1.4 (the Sun has logg = 4.4) to investigate 
the effect of the stellar evolution. In addition, we simulated the stellar winds from a 
3Mq star with log = 2.4 and 1.4. In all the models we use the photospheric magnetic 
field, Brft = 240 G, and the total super-radial expansion factor, / — 240. Note that as 
for the main sequence star (the Sun in this case) these values are between the models 
for fast and slow solar winds (SI06). The stellar wind structures after the quasi-steady 
states are presented in Figure [HI 

The figure shows that hot coronae with T ~ 10 disappears when the star evolves 
to loggr < 2.4, and the density in the atmospheres and winds increase (when measured 
in stellar radii). The increase of the density is owing to the decrease of the surface 
gravity as a result of the expansion of the star; more mass can be lifted up to upper 
regions because of the smaller gravity. The decrease of the temperature can also be 
partly explained by the decrease of the gravity. When the star evolves to the red giant, 
the escape velocity becomes comparable to the sound speed; for example, the sound 
speed of the coronal gas (~ 150 km s~^) exceeds the escape speed of the red giant stars 
with logg ~ 2 at a few stellar radii. Then, the hot corona with T ~ 10 K cannot be 
confined by the gravity and inevitably streams out. 

While the increase of the density is rather continuous with the stellar evolution, 
the temperature rapidly drops from the star with logg = 3.4 to the star with 2.4. 
Thermal instability plays a role in this rapid decrease of the temperature, in addition 
to the gravit y effect explained above. Figure [3 shows the radiative cooling function 
adopted from lLandini fc Monsignori- Fossil ( 1990|) of the solar metallicity gas under the 
optically thin approximation, which we use in the simulation. As shown in the figure, 
the radiative flux decreases for increasing temperature in 10 (K) ~ T ~ 10 (K). 
In this region, the cooling of gas is suppressed when it is heated up; it is thermally 
unstable. Then, the gas is stable in only T ~ 10^ K or T ~ lO*" K in which 
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1 2 5 10 20 

Fig. 6 Evolution of stellar wind from the main sequenee to RGB. From top to bottom, radial 
outflow velocity, Vr (km s~^), temperature, T (K), and density, p(g cm~^), are plotted. The 
solid lines are the results of the 1 Mq stars ; the black, blue, green, and red lines are the 
results of the star with the surface gravity, log 3 = 4.4 (Sun), 3.4, 2.4, and 1.4, respectively. 
The dashed lines are the results of the 3 Mq stars; the black and light-blue lines are the results 
of the stars with the surface gravity, logg = 2.4 and 1.4. 
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T(K) 

Fig. 7 Cooling function of the solar metallicity gas from lLandini fc Monsignori- Fossil lll99Cll V 
The main coolants in the respective temperature ranges are also indicated with arrows. 

thermal conduction also plays a role in the stabilization. This is the main reason why 
the temperature drops from ~ 10 K to ~ 10'^ with stellar evolution. For more detail 
please see S07. 

The properties of the reflection of the Alfven waves are also affected by the stellar 
evolution. Figure [S] compares the reflection of the outgoing Alfven waves between the 
main sequence star (log g — 4.4) and the moderately evolved red giant star (log g = 1.4). 
As clearly illustrated in the bottom panels, the reflection of the outgoing Alfven waves 
is significantly suppressed in the evolved star compared to the main sequence star, only 
30% of the input Alfven waves from the surface is reflected back in the evolved IMq 
star with logg = 1.4, 90% is reflected back in the main sequence star. This is mainly 
because the density slowly decreases in the red giant star and the variation scale of 
the Alfven speed is larger. Then, the outgoing Alfven waves do not sufl^er reflection 
so much and more fraction of the input energy can reach higher altitudes. However, 
because of the larger density the radiative loss is more effective in the evolved stars 
(S07). 

The mass loss rate, A-Kpvrr'^, of the red giant star with logg = 1,4 of 1 Mq star is 
more than 10^ times larger than the mass loss rate of the main sequence star. This is 
much larger than the increase of the stellar surface area (= 1000 t imes); the increase 
of th e mass flux, pvr, itself contributes signiflcantly (S07; see also ISchroder fc Cunt j 
I2OO5I ). 
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Fig. 8 The left panels are the results of the main sequence star with the surface gravity of 
logg = 4.4 and the right panels are the results of the red giant star with log 3 = 1.4 of 1 Mq 
star. The top panels show the wave actions (Sc) of the out-going (solid) and incoming (dashed) 
Alfven waves. The bottom panels show density (solid) and temperature (dashed). The wave 
action {Sc) of the red giant case is reduced to 0.1 times of the original value for comparison. 

The Alfven wave-driven stellar wind plays an important role in terms of the en- 
ergy conversion from magnetic energy to other types of energy. To study the energy 
conversion, a plasma /3 value, which is defined as the ratio of gas pressure to magnetic 
pressure, 

P = Sivp/B^ (9) 

is a useful parameter. Figure |9] shows the /3 values for the different six model^. In the 
region closed to the surface, the /3's decrease on r at first, which is more clearly seen 

^ This figure is modified from the bottom panel of Figure 11 of S07. In S07, the /3 values 
were estimated from the radial magnetic field strength, Br- In this paper, we calculate /3 from 
the total magnetic field strength, the sum of radial and transverse components. 
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Fig. 9 Plasma /3 values as functions of stellar radii. The solid lines are the results of the 
IMq star; the black, blue, green, and red lines are the results of the star with the surface 
gravity, logg = 4.4 (Sun), 3.4, 2.4, and 1.4, respectively. The dashed lines are the results of the 
3 Mq stars; the black and light-blue lines are the results of the stars with the surface gravity, 
logg = 2.4 and 1.4. 



in the unevolved stars. In this region, the atmosphere is mostly static and the density 
decreases rapidly according to an exponential manner. The decrease of the density 
is faster than the decrease of the magnetic energy oc B^, and then, the atmosphere 
becomes magnetically dominated, /? < 1. If the static atmosphere continued to the 
upper altitudes without the dissipation of the Alfven waves, the /3 values would keep 
decreasing. 

In reality, however, the density structure is redistributed by the heating from the 
wave dissipation. In the main sequence (logg = 4.4) and subgiant {\ogg = 3.4) stars, 
the hot coronae are formed, which give larger pressure scale heights ( oc temperature). 
Therefore, the densities decrease more slowly than the magnetic energy, and the j3 
values increase in the coronal regions. In the red giant stars, although hot coronae do 
not form as in the main sequence and subgiant stars, the gas is lifted up directly by 
magnetic pressure associated with Alfven waves in the weak gravity conditions, which 
are gradually connected to the stellar wind regions. Interestingly, the final /3 values 
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Fig. 10 Schematic picture of the local shearing box simulation of an accretion disk. 



are 0.01 — 0.1, being independent from the surface gravities. The redistribution of the 
density through the wave dissipation makes the plasma [3 values stay at the moderate 



4 Turbulent-driven Accretion Disk Winds 

We in vestigated accretion disk winds that are driven by MHD turbulence in accretion 
disks (jSuzuki fc Inutsukalliooil '). The geometrical aspects of disk winds are different 
from those of the stellar winds as easily expected. On the other hand, the driving 
mechanism shows similarities to the solar and red giant winds, which we discuss from 
now. 

To study turbulence driven accretion disk winds, we performed 3D MHD simula- 
tions in a local patch of an accretion disk (Hawley et al. 1995; Matsumoto & Tajima 
1995; Figure I10|l . Here, we neglect the effects of curvature and the simulations are 
carried out in a local Cartesian box that co-rotates with a Keplerian rotating disk. 
Following the convention (Hawley et al. 1995), x, y, and z coordinates corresponds 
to radial, azimuthal, and vertical directions, respectively. While in the x direction we 
do not take into account the radial gradients of the background density, in the z di- 
rection the density stratification owing to the gravity by a central star is included. 
The size of the simulation box is (x, y, z) = {±Q.^H, ±2H, ±411), where a scale height, 
H, is defined from rotation frequency, SI, and sound speed as _ff = \f2cslfl. We use 
(32, 64, 256) mesh points for the [x, y, z) coordinates. The z size of the simulation box 
covers from the sufficiently lower region to upper region of an entire disk so that we 
can simultaneously treat the amplification of magnetic field by MRI and excitation of 
disk winds in a self-consistent manner. This is an advantage of the simulation of the 
accretion disk winds in comparison to the simulations of solar/stellar winds, in which 
we give the conditions at the surface based observationally inferred values without 
considering the generation of the magnetic field in the solar/stellar interiors. 
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Because the co-rotating box is not in an inertial coordinate, inertial forces need to 
be taken into account in the momentum equations : 

^ = -iv(p +^)+ -2^0 XV + sn^x - nlz, (10) 

at p OTT 47rp 

where the third term on the right hand side is Coriolis force, the fourth term corre- 
sponds to tidal expansion which is the sum of centrifugal force and the radial direction 
of the gravity by a central star, and the last term is the vertical component of the grav- 
ity which gives the stratification of the density in the vertical direction. We assume 
isothermal gas instead of solving the energy equation in order to focus of the dynam- 
ics of the disk winds. In the x direction, the shearing boundary condition (Hawley 
et al.l995) is adopted to mimic the differential rotation of a Keplerian rotating disk. 
The simple periodic boundary condition is apphed to the y direction. The outgoing 
condition (SI06) is prescribed to the z direction to properly treat streaming out disk 
winds. We would like to note that our simulation is the first attempt that adopts the 
real outgoing boundary to the disk wind simulation in the local shearing box, whereas 
the simple zero-gr adient boundary was used in previous studies l|Miller fc Stondbood : 



iHirose et al.ll2006h . As the initial condition, we give Keplerian rotation, Vx = —'i/2Qx, 
with small perturbations as seeds for MRI and set up weak vertical magnetic fields 
with /3 = 10® at the midplane. In the vertical direction, we give the hydrostatic density 
structure, p oc exp(— z^/H^). 

In this review paper we mainly focus on the time-averaged structure of disk winds. 
We would like to briefly describe the time evolution and dynamical properties of the 
simulated local accretion disk (see Suzuki & Inutsuka 2009 for more detail). After ~ 3 
rotations, MRI triggers MHD turbulence firstly around z ~ ±(2 — 3)-ff. In the higher 
regions {\z\ > 3-ff), MRI does not set in because the circumstance is magnetically 
dominate from the beginning and stable against MRI. On the other hand, the delay 
of the MRI trigerring in the near-midplane region is an artifact. Since in this region 
the plasma P is too high (gas pressure dominated) at first, we cannot resolve the most 
unstable wavelength of MRI. However, MRI gradually takes place later on. As the 
magnetic field strength increases, the most unstable wavelength can be eventually re- 
solved at the midplane. After 100 rotation, the entire region except the surface regions 
(1 2; I > 3-ff) becomes turbulent and the disk winds are driven by MHD turbulent pres- 
sure from the upper and lower surfaces. At 200 rotations, the magnetic field strength 
almost saturates, balancing the amplification by MRI with the cancellation by mag- 
netic reconnections and the escape with the disk wind^. At this time, the magnetic 
energy is amplified 1000 times of the initial energy of the weak vertical field, and the 
field lines are dominated by the toroidal (y) component (see below). 

An interesting feature is that the disk winds are blow off intermittently. The mass 
fiuxes of the disk winds become strong every 5-10 rotations quasi-periodically. This is 
a consequence of breakups of large scale channel flows at z ~ ±2H (Suzuki & Inutsuka 
2009). 

Figure [TT] presents the disk wind structure averaging over 200 - 400 rotations. The 
variables are averaged on the x — y plane at each z point. The top left panel shows 
that the gas streams out of the upper and lower surfaces. The average outflow velocity 



^ Since the ideal MHD condition is assumed in the simulation, the magnetic reconnections 
take place due to the numerical resistivity determined by the grid scale. Then, the saturation 
level might depend on the resolutions of simulatio ns, whereas this problem is still under debate 
jSano et al.ll2004l : iPessah. Chan, fc Psaitis1l2007f) . 
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Fig. 11 Time-averaged disk structure during t = 200 — 400 rotations. The variables are also 

averaged on x—y plane at each z grid. The top left panel shows Vz/cs (solid), whereas the dotted 
line is the initial condition {vz/cs = 0). The bottom left panel presents the density (solid) in 
comparison with the initial condition (dotted). The top right panel presents magnetic energy, 
B^/47r. The dashed, solid, and dotted lines correspond to x, y, and z components, whereas the 
y component shows both mean (thick) and fluctuation (thin) components. The bottom right 
panel shows the plasma /3 (solid) in comparison with the initial value (dotted). 
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is nearly the sound speed at the upper and lower boundaries. The disk winds are 
dominantly accelerated by magnetic energy (Poyinting flux) of the MHD turbulence, 
which we inspect from now. 

The top right panel shows magnetic energy at the saturated state. The dashed, 
solid, and dotted lines are x, y, and z components. In the y component we are showing 
both mean, (By)'^ and fluctuation, SBy, components. (By)'^ is the simple average on 
X — y planes, {By{z)) — J J dxdyBy{x,y, z)/ (LxLy), and the the fluctuations are 
determined from 5By{z) = J J dxdy(By{x,y, z) — {By(z)))'^ / (LxLy), where Lx{= H) 
and Ly{— 4_ff) are the x and y lengths of the simulation box. As for Bx and Bz the 
fluctuation components greatly dominate the means, and so we simply present Bx and 
Sf. The magnetic energy, which is dominated by the toroidal (y) component as a 
consequence of winding, is amplifled by ~ 1000 times of the initial value (B^ q/Atv = 
10~^) in most of the region {\z\ < 3H). While in the region near the mid-plane {\z\ < 
1.5H), the magnetic held is dominated by fluctuating component (SBy), the coherent 
component {{By)) dominates in the regions near the surfaces {\z\ > 1.5-ff). In the 
surface regions the magnetic pressure is comparable to or larger than the gas pressure 
{j3 ~ 1), and thus, the gas motions cannot control the conflguration of the magnetic 
fields. Therefore, the field lines tend to be straightened by magnetic tension to give 
(5)2 > SB^ there, even if the gas is turbulent. We also note that (B^) is amplified by 
MRI and Parker (1966) instability, whereas (Bz)'^ is strictly conserved. 

The comparison of the final density structure (solid) with the initial hydrostatic 
structure (dotted) in the bottom left panel shows that the mass is loaded up to the 
onset regions of outfiows from z ~ ±2_ff by the amplified magnetic energy. The plasma 
/3 value in the bottom right panel is a good indicator to understand the locations of the 
mass loading. The panel shows that the wind onset regions correspond to /3 « 1, which 
indicates that the winds start to be accelerated when the magnetic pressure exceeds the 
gas pressure. In the further upper regions, \z\ ~ 3H, fi stays almost constant slightly 
below unity. Without disk winds, P decreases with increasing height in the hydrostatic 
structure. However, the mass loading and disk winds inhibit the decrease of p. We 
can interpret that the density structure is redistributed to give the almost constant 
P in the wind regions. This is, in a qualitative sense, similar to the /3 structures in 
the stellar winds (Figure (9]). The detailed features are different mainly because we do 
not solve the energy equation with radiative cooling and conduction in the disk wind 
simulations. 

The disk winds are driven by Foynting fiux; both magnetic pressures and tension 
almost equally contribute. The magnetic tension term, —BzVj_B± (_L denotes x and y) 
shows an interesting distribution, which is shown in Figure [T^ In the surface regions, 
the absolute values decrease with increasing z, because the Foynting flux is converted to 
the kinetic energy of the disk winds; this is natural behavior. However, at \z\ ~ ±1.5-ff, 
the solid line crosses the line. This indicates that the Foynting flux of the tension force 
directs toward both surface and midplane f rom here. This is a resul t of the breakups 
of large scale channel flows at \z\ ~ ±1.5H ( Suzuki fc Inutsukall2009l ). 



The Foynting flux associated with Alfven waves is expressed as the same shape, 
—BzVj_B^, because they are propagating by tension restoring force. The properties of 
the Foynting flux of the magnetic tension shown in Figure [12] are different from the 
Alfven waves discussed in the solar and stellar winds, because dBj_ ^ Bz in the disk 
wind situation. The Foynting ffux here is mainly by large-scale channel flows in mid 
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Fig. 12 Time-averaged Poynting associated magnetic tension, —BzV±B±, flux of magnetic 
tension. 

altitude regions {z ~ 2H) and is associated with a U-shape magnetic field by Parker 
instability in the surface regions (1^1 ~ 3-ff: Suzuki et al.2010). 

5 Summary 

We have introduced the results of our simulation studies on Alfven wave-driven winds. 
In the section[2]we summarized the solar wind simulations presented in SI05 and SI06, 
emphasizing the reflection of Alfven waves. In the standard run, 15% of the initial wave 
energy flux from the photosphere can transmit into the corona, which is sufficient for 
the heating and acceleration of the solar winds. The Alfven waves nonlinearly generate 
compressive waves which eventually steepen and heat up the ambient gas by the shock 
heating. In another way of the explanation, the outgoing Alfven waves dissipate by 
parametric decay to incoming Alfven waves and outgoing compressive waves. This 
process is very efficient in the density stratified atmosphere. 

The structure of the solar winds shows very sensitive dependence on the input wave 
amplitude by the positive feedback mechanism: When increasing the amplitude, the 
coronal density becomes larger and the pressure scale height becomes longer because 
of the larger heating. The larger density enhances the wave dissipation because the 
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wave nonlinearly is larger owing to the smaller Alfven speed (oc l/^fp). The longer 
pressure scale height reduces the wave reflection. Therefore, the heating is enhanced 
when the surface amplitude slightly increases. This shows that the structure of solar 
wind is largely altered on a small change of the input energy through the reflection 
and the nonlinear dissipation of the Alfven waves. 

We introduced recent HINODE observations of Alfvenic oscillations. The obtained 
data bv lFuiimura fc Tsunetal ( 20091) can be explained by the sum of upgoing Alfven waves 



and the comparable amount of the reflected component. They can estimate the leakage 
Poynting flux of the Alfvenic oscillations to the upper location from the phase corre- 
lation of the perturbing magnetic field and the velocity. They analyzed one case and 
found that the sufficient energy to heat the corona and solar wind is going upward. 
These observed properties are quite consistent with what we get in the simulations. 

We discussed the roles of Alfven waves in the context of stellar evolution. Before 
reaching AGB phase when dusts are formed in cool atmosphere, it is expected that, 
similarly to the solar wind, the Alfven waves play a major role in driving the stellar 
winds. By using the surface fiuctuation strength estimated from the convective flux, we 
performed the MHD simulations. When a star evolves to RGB phase, the hot corona 
suddenly disappears and the atmosphere streams out as the cool chromospheric winds 
in which magnetized hot bubbles intermittently float. The main reason of the disappear- 
ance of the hot corona is the gravity effect; the atmosphere flows out before heated 
up owing to the weak gravity confinement. In addition, the thermal instability also 
plays a role in the sudden drop of the temperature. In evolved stars, the Alfven waves 
do not suffer reflection because the density decreases more slowly and the change of 
Alfven speeds become more gradual. This is a part of the reason why the mass loss 
rate jumps up with stellar evolution. 

We also introduced the accretion disk winds by MHD turbulence amplified by MRI. 
The disk winds are driven by Poynting flux when the magnetic energy dominates the 
gas energy. In the disk wind regions the plasma /3 value is kept slightly below unity by 
the redistribution of the density structure, which is qualitatively similar to what we 
got in the simulation of the solar and stellar winds. On the other hands, the disk winds 
are intermittent mainly because of the breakups of large scale channel flows. In the 
simulations, we only model the onset region of the disk winds. In further upper region 
with magnetically dominated condition, Alfven waves might play a role in accerating 
the winds. 
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